Searching for star-planet magnetic interaction in CoRoT 
observations* 
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Abstract Close-in massive planets interact with their 
host stars through tidal and magnetic mechanisms. In 
this paper, we review circumstantial evidence for star- 
planet interaction as revealed by the photospheric mag- 
netic activity in some of the CoRoT planet-hosting 
stars, notably CoRoT-2, 4, and 6. The phenomena are 
discussed in the general framework of activity-induced 
features in stars accompanied by hot Jupiters. The 
theoretical mechanisms proposed to explain the activ- 
ity enhancements possibly related with hot Jupiters are 
also briefly reviewed with an emphasis on the possible 
effects at the photospheric level. The unique advan- 
tages of CoRoT and Kepler observations to test these 
models are pointed out. 

Keywords stars: planetary systems - stars: activity - 
stars: late-type - stars: magnetic fields - stars: general 



1 Introduction 



Massive planets in c lose orbits (semimajo r axis a < 



Kashvap et all boosl Sect. 3.3 



0.15 AU; see, e.g. 
for a motivation for this choice) around late-type 
stars can interact with their hosts through tidal as 
well as magneti c mechanisms (see, e.g., the pioneer- 



ing work by ICuntz et al. I I20OT ). Although the 
study of the tidal interaction is of fundamental im- 
portance to understand th e formation and evolution 
of these systems (see, e.g., iDobbs-Dixon et al. 112004 : 
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Zahn 20081 : Lanza et al. 2011a[) . here we shall focus 
on the interaction involving magnetic fields and the as- 
sociated activity features. The first evidence for such 



an interaction was re ported b v Shkolni k et al. I (j2003l ) 
and later by Shkol nik et al. I ( i2005 , 20081 ). We shall 
review the phenomena observed in some of the hosts 
of the CoRoT planets that may suggest some kind of 
magnetic connection between the stars and their hot 
Jupiters. Since the observations are not conclusive yet, 
further studies are recommended to confirm and clar- 
ify the interaction. We stress the possibilities offered 
by high-precision space-borne photometers, like CoRoT 
and Kepler, for this kind of studies and briefly review 
the mechanisms suggested to explain how a close-in 
massive planet may affect the activity of its host star. 



2 Observations of star-planet magnetic 
interactions 



Shkolnik et al. I (|2003D and IShkolnik et al. I (|2005[ ) ob 



served a sample of twelve hot Jupiter systems reporting 
some evidence for a stellar activity modulated with the 
orbital period of the planet rather than the star rotation 
period. The cases of HD 179949 and v And are those 
displaying the clearest effects in observations of modu- 
lated Ca II H&K chromospheric emission. Their chro- 
mospheric hot spots rotating in phase with the planet 
irradiate powers up to ~ 10^^ W and are not located 
at the subplanetary longitude, but lead the planet by 
« 80° in the case of HD 179949 and w 160° in the case 
of V And. In both the cases, the hot spots seem to 
have a lifetime of ^ 300 — 400 days and are detected 
only in ^ 50 — 60 pe rcent of the seasons. Therefore, 



Shkolnik et al. 



()2008h suggested that the star-planet 
magnetic interaction (hereafter SPMI) undergoes tran- 
sitions fror n. "on" to "off' states and vi ceversa. Obser- 
vations by Poppenhaeger et al. I ( 201l[ ) did not detect 
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any hot spot modulated in phase with the planet in 
V And during their campaign in 2009. 

The signatures of S PMI were searched als o in the 
coronal X-ray emission. IKashvap et al. 1 (I2OO8I) claimed 
that stars with hot Jupiters (a < 0.15 AU) have a 
higher X-ray flux than stars with distant (a > 1.5 AU) 
planets, with an enhancement by a factor of w 2. 
However, Poppenhaeger et al. ( 2010l ) cast doubts on 
this result and found no evidence of an X-ray emis- 
sion enhancement by comparing volume-limited sam- 
ples of stars with and without hot Jupite rs. They ex- 



plain the result of Kashvap et al. ( 20081) 



and HD192263. IWalker et al. I pOOSi ) observed t Boo 
for five years both from the ground by monitoring the 
flux in the core of the chromospheric Ca II H&K lines 
and from the space with the MOST (Microvariability 
and Oscillation of Stars) microsatellite. A persistent 
active region was found on the star at a longitude lead- 
ing the subplanetary longitude by 80° with a variable 
brightness and a lifetime of at least w 1400 days. Since 
the rotation of t Boo is almost synchronized with the 
orbital motion of its planet, it was the persistence of 
the active region that pointed towards a SPMI effe ct. 



as a conse- 



quence of the sele c tion e ffects affecting their s amples. 
Claims by IScharf (I2OIOI) and iHartman I (|2010l) about 
a relationship between stellar X-ray luminosity and the 
mass or the surface gravity of the planet, respectively, 
were also refuted as consequences of subtle selection ef- 
fects affecting the samples used by th ose investigators 
(see IPoppenhaeger &: Schmitt II2OIII and the contri- 
bution by Katja Poppenhaeger to the 2nd CoRoT Sym- 
posium). However, an iiiteresti ng observation was re- 
ported bv lPillitteri et al. I |2010l ) who found a softening 
of the X-ray spectrum and a strong flare close to the sec- 
ondary eclipse of the hot Jupiter orbiting HD 189733. 
Although this could simply be a coincidence, a coronal 
flaring activity modulated with the orbit of the planet 
cannot be excluded. 

Interesting results came from observations in the 
ultraviolet, performed mainly wi th the Hubble Space 



In the case of HD 192263, ISantos et al. I (|2003l ) re- 
ported evidence of starspots rotating with the planet 
orbital period for several rotations and lagging the sub- 
planetary longitude by « 90°. This may produce a con- 
fusion between the radial velocity jitter induced by ac- 
tivity and the actual reflex motion of the star produced 
by the planet since both have the same period, but the 
former lacks the long-term phase coherence that char- 
acterises the latter. Further evidence of photospheric 
features possibly associated with SPMI has been ex- 
tracted from CoRoT photometry, as we discuss in the 
next section. 

For completeness, we mention also the possibility of 
detecting SPMI in the radio domain. The method is 
promising because it can allow us to directly mea- 
sure the planetary magnetic field (cf. end of Sect. 14. 1) ). 
although no confirmed radio de tections of exoplan 



Telescope. The first studie^.g., I Vidal-Madjar et al. | iLecavelier Pes Etangs et al.1l20lll) 



2003; Ehrenreich et al. 20081) found an extended hy- 
drogen envelope around the transiting planet HD 209458b| 
likely to be the result of the evaporation of the plan- 
etary atmosphere under the act ion of the stellar ir- 



ets have been reported to date dHess fc Zarka 2011 



radiation. More recent studies ( Fossati et al. 2010l) 



detected an excess absorpti on in the nea r UV during 
transit ingress in WASP-12. iLlama et al. I (|201ir) mod- 
elled it as resulting from the absorption of the stellar 
coronal plasma compressed at a bow shock formed at 
the boundary with the planetary magnetosphere. Al- 
though the present observations are still preliminary, 
UV transit light curves can offer a unique opportunity 
to d etect such bow s h ocks, as theoretically discussed by. 



e.g 



Vidotto et al. 



(|2011l ). The shock characteristics 
depend on the geometry and intensity of the planetary 
magnetic field as well as on the properties of the stellar 
coronal plasma and the stellar wind. Therefore, they 
can be used to derive information on the planetary mag- 
netic fields given that the stellar coronal field can be 
modelled through the extrapolation of the photospheric 
field measured by spectropolarirnetric techniques (e.g. 



3 Photospheric SPMI from CoRoT 
observations 

We introduce three cases, i.e., those of CoRoT-2, 
CoRoT-4, and CoRoT-6. In the case of CoRoT-2, the 
mean rotation period of the star, P^ot = 4.522 days, 
is longer than the orbital period of the planet, Porb — 
1.743 days. On the other hand, CoRoT-4 is synchro- 
nized with a mean rotation period of ^ 8.8 days, 
very close to the orbital period of 9.202 days, while 
CoRoT-6 is an example with an orbital period (Porb = 
8.889 days) longer than the mean rotation period of 
the star (Piot = 6.35 days). We refer to the mean ro- 
tation period of the stars because they show surface 
differential rotation with a mini mum amplitude of 3 
perce nt in the case of CoRoT-2 (ISilva-Valio fc Lanza 



sp' 

Moutou et al~l 12007! : liares et al. Iboioh 



20111 1 , ~ 6 percent for CoRoT-4 jlLanza et al. 
and ~ 12 percent for CoRoT-6 ([Lanza et al. 



2009 



2011 



SPMI signatures in the photospheres of planet- 
hosting stars have been suggested in the case of t Bootis 



These cases are interesting also in view of the different 
methods applied to extract the signal associated with 
SPMI according to the ratio Prot/-Porb- 
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3.1 CoRoT-2 



The light curve of CoRoT-2 has been exten sively mod- 



Lanza et al. 



ja) ap- 



elled by several researchers, 
plied a maximum entropy spot model to the out-of- 
transit light curve to derive the distribution of the cov- 
ering factor of the active regions vs. longitude along 
successive rotations. A comparison of this approach 
with other spot inodelli n g methods is provide d by. 
e^ 



Mosser et all (l2009l) . ]Frohhch et al7l (|2009l ). and 



Huber et al. 



while a detailed test with solar 



observations is presented by iLanza et al. I (|20Q7I ). The 



modelling technique exploits the modulation of the vis- 
ibility of the active regions by the rotation of the star, 
so that a good phase coverage of each rotation is needed 
to derive a stable and reliable map. The advantage of 
CoRoT is the almost perfect time sampling (the duty 
cycle is virtually 100 percent for spot modelling), but 
the duration of the stellar rotation sets a minimum 
timescale for an accurate mapping that is ^ 3 days in 
our case, i.e., about 65 percent of a rotation. Morover, 
since the active regions evolve rapidly, it is not possible 
to obtain an adequate best fit for time intervals longer 
than 3 — 4 day s. In conclusion, the time resolution of 



the mapping of lLanza et al. I (|2009a|) is longer than the 



orbital period of the planet making impossible a direct 
search for SPMI features in those maps. 

The distribution of the active regions versus longi- 
tude and time on the surface of CoRoT-2 is plotted in 
Fig. [1] where we see two persistent active longitudes in 
which individual active regions appear, grow, and de- 
cay with lifet imes up to ~ 50 — 60 days. The model by 
Lanza et al. ] (j2009a,) assumes that each active region 
consists of dark spots and bright solar-like faculae with 
a facular-to-spotted area ratio indicated by Q. The 
models in Fig. [T]have been computed for Q = 0.0, i.e., 
without faculae. The variation of the total spotted area 
vs. time is plotted in Fig. [2] and shows oscillations with 
a period of 29.6 ±4.0 days, both in the model assuming 
only dark spots as well as in that with a facular com- 
ponent having Q — 1.5. Note that for the Sun Q — 9, 
while for more active stars, such as CoRoT-2, lower 



values of Q are expected (jLanza et al. 1 120071 ). 



Two possible interpretations have been suggested for 
the oscillations of the total spotted area. The first 
considers them analogous to the oscillations of the to- 
tal sunspot area observed close to some of the max- 



ima of the eleven- year cycles i n the Sun (j Oliver et al. 



1998t Zaaarashvili et al. 2010) and displa ying a. peri- 
odic ity of 150 — 160 days. A ccording to iLou (I2OOOI) 



and IZagarashvili et al. 



(l2010l ). Rossby-type magneto- 
hydrodynamic waves, excited in the subphotospheric 
layers or at the interface between the convection zone 




1 CO 20'0 
Longitude (deg) 



Fig. 1 Maximum entropy spot model of CoRoT-2. The 
isocontours of the spot filling factor are reported vs. lon- 
gitude and time. Different colours refer to different values 
of the filling factor with dark blue indicating the minimum 
(virtually no spots) and orange/yellow the maximum. The 
longitude is measured in a reference frame rotating with the 
mean stellar rotation period of 4.5221 days and increases in 
the direction of the stellar rotation a nd the orbital motio n 
of the planet which is prograde (cf. iBouchv et al. |[2008l ). 
Note that the longitude is repeated beyond the 0° — 360° 
interval to he l p follow ing the migration of the spots. Credit: 
iLanza et al. I (|2009al ). reproduced with permission, © ESO. 
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Fig. 2 The total spotted area of CoRoT-2 vs. the time. 
The solid line refers to models with dark spots only (Q = 0), 
while the dotted line indicates the spot area in the mod- 
els with a facular component with Q = 1.5. Credit: 
iLanza et al. I ([20093), reproduced with permission, © ESO. 
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and the radiative interior of the Sun, could aceount for 
such a short-term periodicity by modulating the emer- 
gence of magnetic flux. Since the period of the wave is 
proportional to the rotation period of the star, we ex- 
pect a period of about one month in the case of CoRoT- 
2 that is a G-type star rotating five times faster than 
the Sun. 

An alternative explanation for the modulation of the 
spotted area in CoRoT-2 considers it as a signature of 
SPMI in the photosphere of the star because of the 
close coincidence of the periodicity with ten synodic 
periods of the planet with respect to the mean stellar 
rotation period. The synodic period Psyn is the time 
interval between two successive passages of the same 
meridian across the subplanetary longitude in a rotat- 
ing star. It is equal to 2.89 days in the case of CoRoT-2, 
assuming a mean rotation period Piot — 4.5221 days be- 
cause it is given by: p-^\ = \P;'^^ - P~^}^\. According 



to this conjecture, the passage of the planet over an ac- 
tive region may trigger the emergence of magnetic flux 
tubes when their intensity is already close to t he thresh- 



old for the onset of the buoyancy instability (jAcheson 



19781 I979I) . Ten passages are to elapse before the mag- 
netic field intensity reaches again the threshold needed 
for the planetary-induced perturbation to be effective. 
In Sect, m a possible physical mechanism to account 
for such a perturbation will be introduced, based on 
the effects on the stellar dynamo induced by the re- 
connection in the stellar coronal field triggered by the 
planetary motion. 

Spots occulted by the planet during transit can be 
detected through the characteristic light bumps pro- 
duced along the transit light curve when the planet's 
disc passes over them. The average spotted area 
of the occulted spots plotted vs. longitude shows 
six relative maxima, one of which is associated with 
the subplanetary longitude within ±3° (cf. Fig. 4 of 

where the subpla netary 



Silva-Valio fc Lanza 



2011 



longitude is zero). Silva-Valio h Lanza! ( 2011 ) inves 



tigated the variation of the total area of the occulted 
spots and found a periodicity of 17.7 ± 2.3 days. Al- 
though very close to ten orbital periods of the planet, it 
is unlikely to be an effect of the aliasing due to the pe- 
riodic time sampling of the observations. In that case, 
we should observe other periodicities at, say, 5 or 15 
orbital periods which is n ot the case (see Fig. 10 in 



Silva-Valio &: Lanza 1120111) . Moreover, computing the 



periodograms of 10 000 random permutations of the 
time series, a periodicity of 17.7 days occurs only in 
~ 3 percent of the cases. 

The periodicity corresponds to ~ 6 synodic periods 
when the rotation period of the mean latitude of the 
strip occulted by the planet is considered, i.e., 4.48 days 



(Sil va-Valio h, Lanza 1120111 ). The occulted spots are a 
subset of all the spots present on the star and are lo- 
cated at a latitude between ~ 5° and ~ 20°. On the 
other hand, the modulation of the light curve outside 
transit is produced by starspots at all latitudes. We 
expect that spots can be formed in a latitude range 
more extended than in the Sun because CoRoT-2 is at 
least 20 times more active than our star as measured 
by the amplitude of the optical fl ux modulation pro - 



Strassmeier 20091) 



duced by the active regions (e.g. 
Therefore, the periodicity of the starspot area of the 
band occulted during transits is diluted when we con- 
sider the starspot area derived from the out-of-transit 
light curve. As a matter of fact, a peak correspond- 
ing to the second harmonic of the 17.7 day period is 
present in the periodogram of the out-of-transit spot- 
ted area although its height is remarkably lower than 
that of the main peak at 29 days because of the e ffects 



of the dilution (cf. Fig. 7 of iLanza et al. Il2009al) . 



A possible interpretation of the 17.7 day periodicity 
that corresponds to ^ 6 synodic periods, is that the pas- 
sage of the planet over the low-latitude active regions 
in the occulted band induces the same phenomenon as 
observed in the variation of the full-disc area. How- 
ever, its cadence is shorter because the magnetic field 
is amplified more rapidly at low latitudes. This inter- 
pretation assumes that there are two latitudes where 
the conjectured mechanism is operating, one in the oc- 
culted band, the other in the latitude range not occulted 
by the planet. The localization of the triggering in two 
latitude bands is required by the absence of interme- 
diate periodicities in the area modulation. Moreover, 
these periodicities appear only when the total spotted 
area, i.e., integrated over the longitude, is considered. 
Therefore, the modulation is not a property of some 
specific active region, but of all the regions in a given 
latitude interval. 



Finally, iPagano et al. 



(|2009f ) considered the vari- 
ance of the flux of CoRoT-2 in the CoRoT white pass- 
band vs. the orbital phase. The variance was calculated 
using the individual 32 s exposures binned vs. the or- 
bital phase. Considering the first 75 days of the light 
curve when the jitter effects were minimal and a bin 
size of 0.05 in phase, they found that the variance of 
the flux varied regularly with the orbital phase reach- 
ing a maximum just immediately before the transit and 
a minimum around phase 0.4. This phenomenon could 
be related to an enhancement of the activity in CoRoT- 
2 triggered by the passage of the planet, e.g., in the 
form of small white-light flares, maybe associated with 
the spots generally found at the subplanetary longitude 
during transits (see above). 
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3.2 CoRoT-4 

In the case of CoRoT-4, the mean rotation period of 
the star is close t o the orbital p e riod of the planet. The 
spot modelling of Lanza et al. ( 2009bl ) shows a persis- 
tent active region located ±50° from the subplanetary 
longitude, when a longitude reference frame rotating 
with the orbital period of the planet is assumed (see 
Fig. [3]). Note that the facular-to-spotted area ratio Q 
adopted for the modelling affects somewhat the distri- 
bution of the spotted area vs. longitude. This happens 
because faculae are brighter when they are closer to the 
limb and thus they affect the distribution of the dark 
spots to reproduce the o bserved light modulation (cf., 



e.g., iLanza et al. Il2007t) . Nevertheless, the presence of 



spots within ±40° from the subplanetary longitude is 
confirmed by the modelling of the light bumps associ- 
ated with the spots occulted during the transits (Silva- 
Valio & Lanza 2011, in prep.). 

The star shows a differential rotation of at least 
6 percent which is not unusual for an F-type star. 
In spite of that, the active longitude associated with 
the planet has persisted for at least ^ 60 days, i.e., the 
duration of the CoRoT observations. Unfortunately, 
this time span is too short to exclude that the coinci- 
dence happened by chance, but the phen omenology is 



reminiscent of what IWalker et al. I (|2008l ) observed in 



r Boo, another F-type star in almost synchronous ro- 
tation with the orbit of its hot Jupiter. In that case, 
in spite o f a remarkable diffe rential rotation (~ 18 per- 



cent; see 



Catala et al.ll2007() , a persistent active region 
was found leading the planet by « 80° with a lifetime of 



at least « 1400 days. iMosser et all (|2009„ ') found that 



photospheric spots have typical lifetimes not exceeding 
one week in F-type stars which makes such a long per- 
sistence remarkable. The possibility that the feature is 
associated with the tidal interaction is ruled out by the 
lack of a similar feature on the opposite hemisphere of 
the star, i.e., associated with the other tidal bulge. 

3.3 CoRoT-6 

A maximum entr opy spot mod e Uing of CoRoT-6 has 
been presented by Lanza et al. ( 2011bl ) and shows an 
interesting phenomenology. Now the rotation period 
(-Plot = 6.35 days) is shorter than the orbital period 
of the planet (Porb = 8.886 days), giving us the possi- 
bility to resolve the spot evolution during the orbital 
as well as the synodic periods (Psyn = 22.25 days). 
Relative maxima of the spotted area occur in about 
ten active regions when a longitude lagging the planet 
by ~ 200° passes by the regions themselves, as shown 
in Fig. m The probability of a chance association has 



been estimated by Lanza et al. ( 2011bl ) to be lower 
than 1 percent. Again, the passage of the planet over 
the stellar active regions seems to be associated with 
a triggering of new magnetic flux emergence, although 
there is an average time lag of ^ 3.5 days between the 
passage of the planet and the maximum spotted area 
in the active regions. The total spotted de- 
rived from the out-of-transit light curve, does not show 
any obvious periodicity, possibly because the synodic 
period is comparable or longer than the lifetime of the 
active regions, so the planet-induced modulation is hid- 
den by the larger variations associated with the intrin- 
sic growth and decay of the active regions. On the 
other hand, a preliminary analysis of the variation of 
the total spotted area occulted during transits shows a 
modulation with a period of 22.5 — 24 d, close to the 
synodic period of the planet (Silva-Valio & Lanza 2011, 
in prep.). 

3.4 Comparison with some preliminary Kepler results 

The preliminary results coming from the analysis 
of the light curves of the Kepler space telescope 



( Borucki et al. I 120091 ) show two interesting cases of 



possible star-planet int e ractio n in the photosphere 



Sanchis-Oieda fc Winn I (|201l[ ) analysed the transit 
light curves of HAT-P-ll, a K4V star, confirming that 
the system is strongly misaligned as already known 
from the observations of the Rossiter-McLaughlin ef- 
fect. The interesting feature is that the spot occul- 
tations by the planet occur only at two fixed phases 
during the 26 observed transits. This is explained by 
assuming that the planet transits over the starspots lo- 
cated in two low-latitude activity belts on a star viewed 
almost equator-on, similar to what we observe in the 
Sun where sunspots appear in two active latitudes sym- 
metric with respect to the equator that migrate towards 
the equator as the cycle progresses. Another possibility 
is that the star is viewed almost pole-on, in which case 
the transit chord cross es twice the same high-latitud e 
spot belt (cf. Fig. 7 in ISanchis-Oieda fc Winnlboill) . 
The intriguing feature is that starspots are observed at 
the first intersection between the transit chord and the 
active belts in 13 out of 26 transits and at the second 
intersection in 13 out of 26 transits. In other words, 
in the 50 percent of the transits we observe a starspot 
not too far from the subplanetary longitude. However, 
since the star is about five times more active than the 
Sun at the maximum of the eleven-year cycle, as mea- 
sured by the amplitude of the optical light curve, this 
close association between the planet and the occulted 
spots could have occurred simply by chance. 
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Fig. 3 Averaged spot filling factor vs. longitude in 
CoRoT-4. The longitude is measured in a reference frame 
rotating with the orbital period of the planet and increases 
in the direction of the stellar rotation and the planetary 
motion, assumed to be prograde. The top panel shows the 
spot model with a facular component in the active regions 
(Q = 4.5), while the lower panel shows a model assum- 
ing only dark spots (Q — 0). The vertical dashed lines 
mark the longitudes 0° and 360° beyond which the plot is 
repeated, while the verti cal dotted line marks the subplan- 
etary longitude. Credit: iLanza et al. I (|2009bl ). reproduced 
with permission, © ESO. 



A more impressive case is that of Kepler- 17, a G- 
type star quite similar to the Sun, except for its rota- 
tion period of 11.9 d, and accompa nied by a hot Jupite r 
with an orbital period of 1.486 d (jPesert et al.ll201ll) . 
The rotation period of the spots occulted during the 
~ 180 observed transits is just eight times the planet 
orbital period producing a particular stroboscopic effect 
in their visibility. Five occulted active regions show a 
lifetime of at least ~ 100 d and reach their maximum 
filling factor when they cross the subplanetary longi- 
tude corres ponding to phase in the right hand plots 
in Fig. 11 of lDesert et al.l (|201l[ ). This remarkable asso- 
ciation between the planet and the occulted starspots 
is confirmed by the average transit light curve which 
shows a clear bum p centred at phase (cf. Fig. 2 in 
Desert et al.ll201l[ ). In view of the large number of ob- 
served transits covering a time interval of 16.7 months, 
such an association is highly unlikely to have occurred 
by chance. Moreover, the phenomenology is remarkably 
similar to that of CoRoT-6, in spite of the differences 
in the stellar and planetary parameters. 



4 Theoretical models 



4.1 Chromospheric and coronal features of SPMI 




1 00 200 
Longitude (deg) 



Fig. 4 Maximum entropy spot model of CoRoT-6 with 
Q = 1.5. The isocontours indicate different values of the 
starspot filling factor with the same color coding of Fig. [T] 
The adopted longitude reference frame is rotating with the 
mean rotation period of the star (Prot ~ 6.35 days). The 
longitude increases in the direction of the stellar rotation 
and the orbital motion of the planet (which is prograde) , and 
is repeated beyond the 0° — 360° interval to help following 
the migration of the spots. The white crosses mark the rel- 
ative maxima of the spot filling factor in the active regions, 
while the white straight lines indicate a lon gitude lagging 
the sub planetary longitude by 200°. Credit: iLanza et al. I 
reproduced with permission, © ESO. 



Several theoretical models have been proposed to ac- 
count for specific features of SPMI as observed in the 
chromosphere and possibly in the corona, but a compre- 
hensive theory is still lacking. To explain the presence 
of chromospheric hot spots rotating in phase with the 
planet, the first idea was that of considering a scaled 
version of the unipolar induction model proposed for 
the Jupiter-Io system. Observations in the UV has re- 
vealed two bright spots in the Northern and the South- 
ern hemispheres of Jupiter, respectively, located at the 
footpoints of the loops connecting the magnetic poles 
of Jupiter with lo. A flux of Alfven waves is excited 
by the motion of lo across the magnetic field lines of 
Jupiter's large-scale dipole field and the energy of the 
waves is conducted down t o the loop foo tpoints produc- 
ing the emission (see, e.g., IZarka 1120071 and references 
therein). This mechanism can work even if lo has no 
intrinsic magnetic field because it is sufficient its orbital 
motion across the Jupiter's field lines and some surface 
electric conductivity to excite the Alfven waves. 

A limitation of the model in the stellar case is the 
low power that can be conveyed onto the chromosphere, 
estimated of the order of 10^^ — 10^'' W. The observed 
phase lag between the planet and the hot spot is also 
difficult to explain wit h a potential dipole magnetic 
field like that of Jupiter. iMcIvor et al.l (j2006il suggested 
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that the axis of the stellar dipole is tilted with respect 
to the orbital angular momentum of the planet, but 
even with this hypothesis it is not possible to account 
for a phase lag of ~ 160° as observed in the case of 



V And . On the other hand, the model bv lPreusse et al. 



( 20061). recently confir med by the numerical simulations 



of 



Kopp et al. 



(|201lh . can explain a large phase lag. It 
is based on packets of Alfven waves that are excited by 
the orbital motion of the planet, even in the absence of 
an intrinsic planetary field, and then propagate towards 
the star along characteristics making some angle with 
the magnetic field lines. They can reach the surface 
provided that the velocity of the stellar wind at the dis- 
tance of the planet is subalfvenic, that is the case when 
the planet is sufficiently close to the star. The main 
limitation is the low energy flux at the stellar surface 
that makes almost impossible to account for an emitted 
power as high as 10^° — 10^^ W. Note that the same 
order of magnitude is required to expl ain the enhance- 



ment of the X-ray emission claimed by iKashvap et al. 
( 20081) . 



To explain the emitted power, iLanza I (|2009l ) con- 
sidered the possibility that the interaction between the 
planetary magnetic field and the stellar coronal field 
may trigger the release of the energy already stored 
in the coronal loops. The reconnection process con- 
tinuously produced by the planetary field in the outer 
corona tends to reduce the magnetic energy of the coro- 
nal field. Since the minimum energy state for a given 
total magnetic helicity is a linear force-free field, the ge- 
ometry of the field lines can be described analytically 
given the boundary conditions at the stellar surface. 
An interesting property of linear force- free fields is that 
of being twisted, thus providing a natural explanation 
for the phase lag between the planet and the chromo- 
spheric hot spot located close to the footpoints of the 
magnetic field lines. Th is property has been explored in 
detail by Lanza ( 20081) . In combination with t he trig- 



gering of magnetic energy release suggested by iLanza 



(|2009f l. brce-free fields could account for the emitted 
power of the hot spots as well as for some coronal 
emission enhancement. Depending on the total mag- 
netic helicity and the surface boundary conditions, the 
magnetic field lines intersecting the planetary magne- 
tosphere may close up before reaching the surface of 
the star which may imply that no magnetic connec- 
tion exists betwee n the planet and the chromosphere 
(see iLanza 1120091 ). Such configurations would give rise 



to an "ofF' state of t he SPMI in the terminology of 
IShkolnik et al. ( 2008 ). When the boundary conditions 
and/or the total helicity change, a field configuration 
that connects the planet with the star can be resumed 
leading to a switching on of the SPMI features in the 



chromosphere. Another interesting property of linear 
force-free fields is the predominance of closed field lines 
that reduce the loss of angular momentum in the stellar 
wind and the consequent braking of the stellar rotation. 
This may have important consequences for the estimate 
of the age of the stars accompanied b y hot Jupite rs 
through the method of gyrochronology ( Lanza 2010l) . 



Full magnetohydrodynamic simulations, such as 



those of llpet al. I I2OO4I ) ICohen et al.l (l2009al) ICohen et aCk 



( 2009bl). 



Cohen et al. 



and Cohen et al. 



(2010D, ICohen et al. I (l2011a 



l2011bl ). confirm several results of 

the a nalytical force- free models by iLanza I (I2OO8L I2QO9I 
2OIOI) and show how a stellar corona tends to be con- 



fined and heated by the interaction with the magne- 
tosphere of a hot Jupiter. Hot coronal regions can 
account for an increased emission in the X-rays and 
the reduced efficiency of the acceleration of the stel- 
lar wind. Although the magnetohydrodynamic regimes 
accessible to numerical simulations are many orders of 
magnitude far from real systems, these studies provide 
a wealth of information on what is to be expected in 
real star-planet interaction and can be specialized to 



the ca se of particular systems (cf., e.g., ICohen et al. 



2011aD 



Another approach to the star-planet interaction has 
considered individual coronal loops reconnecting with 
the planetary magnetosphere. According to the phase 
of the stellar activity cycle, the height of the tops of 
the loops can vary leading in some cases to a strong 
interaction, modulated with the orbital period of the 
planet. On the other hand, during phases when the 
loops are not so tall the interaction with the plan- 
etary field is reduced and we pred ominantly see the 



modu lation by the stellar rotation (jCranmer fc Saar 



20071) . Thus this mechanism can naturally account 



for the on /off t ransit i ons i n the SPMI reported by 
Shkolnik et al~l |2005l I2OO8I ). Other studies have ad- 



dressed the transport of energy from the reconnection 
site, presumably in the corona, to the chromosphere 
through beams of energetic particles that produce a 
localize d heating when they impact onto the chromo- 
sphere |Gu & Suzuki 112009!) . 

Finally, we should mention the models developed 
to predict the radio emission from the systems host- 
ing hot Jupiters. Since those planets are much closer 
than Jupiter to the Sun, the stellar wind speed is likely 
to be subalfvenic. However, a remarkable interaction 
between the wind and the planetary magnetosphere is 
predicted, and the emitted radio fluxes should be of the 
order of 10^ those of Jupiter. The peak of the emitted 
spectrum should fall at frequencies of a few or several 
tens of MHz depending on the strength of the plan- 
etary magnetic field (see, e.g., Stevens 2005t Zarka 
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2007; Jardinc & Collier Cameron 2008: Lanza 2009i 



Hess fc Zarka and references therein). A signif- 

icant flux at microwave frequencies is expected if high- 
energy particles spiral along the kC magnetic fields of 
the spots in the stellar photosphere, but to attribute 
the emission to SPMI it is necessary to detect its mod- 
ulation with the orbital or the synodic period. On the 
other hand, radio emission in the MHz range is charac- 
teristic of the planetary fields and its detection would 
provide information on the strength of the planetary 
magnetic fields. 

4.2 Accounting for the energy budget 

Most of the models considered above dealt with analo- 
gies or differences between the solar-planet and the star- 
planet interactions, or focussed on specific issues such 
as the phase lag between the chromospheric hot spot 
and the planet. Nevertheless, only some of them con- 
sidered the problem of dissipating a power of the order 
of 10^° — lO'^^ W to account for the energy irradiated b y 



the chromospheric hot spot s of IShkolnik et al. I (|2005[ ). 



In their MHD simulations, lip et al. I (|2004l ) were able 
to obtain powers of the order of 3 x 10^^ W adopt- 
ing a magnetic field of 0.1 G and a size of the inter- 
action area of 5 Jupiter radii at the magnetopause be- 
tween the planetary field and the stellar coronal field 
with an orbital velocity of the planet of 280 km s~^. 
However, typical surface fields in HD 179949 are of the 
order of 10 G at most and even in the case of the re- 
markabl y active host HP 18 9733 they do not exceed 
- 40 G (jMoutou et"al~ll2007t) . Since the magnetic field 
strength decreases at least as {r/R)~^, where r is the 
distance from the centre of the star and R its radius, we 
expect field strengths B ~ 0.02 — 0.08 G at a distance of 
8R, typical of most of the systems with hints of SPMI. 
Since the dissipated power is proportional to B^, the 
field strength is too low to account for the energy radi- 
ated in the chromospheric hot spots. A power at least 
t wo or three ord ers of magnitude larger is required. 
Lanza ( 2009f ) proposed that the energy needed to 



explain the hot spot comes from the whole magnetic 
structure connecting the planet with the star and is 
released at a closer distance than the planet orbital ra- 
dius. His considerations were based on models bor- 
rowed from solar physics to explain the budget of solar 
fiares and considered the fundamental role played by 
magnetic helicity in the evolution of the coronal field. 
Here, we introduce some intuitive ideas for such a kind 
of models by means of a simple cartoon (see Fig. [S]). 
We consider a late-type star with a hot Jupiter whose 
orbital motion is not synchronous with stellar rotation. 
At the beginning of the process, a long coronal loop 



reconnects with the planetary magnetosphere releasing 
a modest amount of energy, with a characteristic dissi- 
pated power of the order of 10^^ W for a time scale of 
the order of L/va, where L is the size of the interaction 
region and va the Alfven velocity. Once the reconnec- 
tion process is over, the relative orbital motion of the 
planet with respect to the loop footpoints produces a 
bending of the loop on a timescale of the order of the 
day. The stress of the field accumulates energy into the 
loop and the field is no longer potential as indicated 
by the curvature of the field lines in the outer corona. 
This energy can be released closer to the star when the 
curved field lines interact with shorter loops, as in the 
final sketch of the cartoon. Since the interaction site is 
closer to the star, the magnetic field intensity is remark- 
ably larger there, while the relative velocity field, that 
scales roughly as {r/R), is still of the order of several 
tens of km s~^. Considering a field strength B = 1 G 
at the interaction site, a typical size of the interacting 
region oi L = R, and a relative velocity v — 30 km s^^, 
we have a released power of: 



Pdiss ^ j-B'L'v ~ 3.5 X 10"% W, 



(1) 



where < 7 < 1 is a factor depending on the relative 
angle between the reconnecting magnetic field lines at 
the interaction site and fi is th e magnetic pe rmeability 
of the plasma (cf. Eq. (8) of iLanza 1120091 ). For op- 
positely directed fields, we get the maximum efficiency 
with 7 ~ 1 and we have a power sufficient to explain 
the energy radiated by chromospheric hot spots. Note, 
however, that ordinary stellar flares release a compara- 
ble power when the relative velocity between magnetic 
loops is of the same order. The effect of the planet is 
that of inducing a sequence of fiares when an intercon- 
necting loop is stressed again by the relative motion of 
the planet. Since most of the energy is released close 
to the star, the chromospheric heating does not require 
any ad hoc process to convey the energy from the plan- 
etary orbit to the chromosphere. The only requirement 
is the presence of an interconnecting loop whose field 
lines are stressed by the relative motion of the planet 
without producing an instability that destroys the loop 
itself before it have had time to accumulate enough en- 
ergy. An important factor is the transverse size of the 
interconnecting loop, that we assume of the order of the 
radius of the star. It depends also on the intensity of 
the magnetic field of the planet because a stronger field 
pro duces a more extended magnetosphere (cf. Eq. (7) 



Lanza 20091 ) that increases the probability of re- 



connection with stellar loops and thus the transverse 
size of the interconnecting loop. The energy extracted 
from the orbital motion and stored in the magnetic in- 
terconnecting loop should in principle lead to a slow 
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orbital decay of the planet. Since the total mechani- 
cal energy of the orbital motion is GMsmp/2a, where 
G is the gravitation constant, Ms the mass of the star, 
and rUp that of the planet, a typical orbital energy of 
Iq36 _ |q37 j jg obtained for a hot Jupiter. With the 
dissipated powers considered above, this implies an or- 
bital decay timescale of at least several Gyr. There- 
fore, the tidal interaction is likely to be more effective 
for the evolution of t he orbit of the planet (cf., e.g., 
Goupil fc ZahTlboOSl) . 

This promising model for the magnetic star-planet 
interaction, here sketched with simple order-of-magni- 
tude and qualitative considerations, will be presented 
in more detail in a forthcoming study (Lanza 2011, in 
prep.). Note that in systems close to synchronizations, 
such as T Boo or CoRoT-4, the relative velocity v is 
expected to be remarkably smaller than in the case 
discussed above. However, surface differential rotation 
may still provide sufficient relative velocity for the onset 
of the process, although with a reduced efficiency. 

4.3 A conjecture on the perturbation of the stellar 
dynamo by the planet 

The photospheric phenomena presented in the case 
of the CoRoT and Kepler targets in Sect. |3] or in 
Walker et al. I ( 2008^ for the case of r Boo cannot be ex- 



plained with the models discussed so far because they 
consider only the interaction and the associated phe- 
nomena in the outer atmosphere. Nevertheless, as we 
shall see in this section, there could be a connection 
between the perturbation produced by the planet on 
the magnetic field of the corona and the operation of 
the stellar dynamo responsible for the amplification and 
modulation of the stellar field appearing in the photo- 
sphere. 

The emergence of new magnetic fiux in the form 
of spots or faculae requires some mechanism to trig- 
ger the buoyancy instability of the magnetic field when 
the planet passes over an active region. The first 
conjecture is that the tidal deformation of the stel- 
lar field may promote its instability if the field is a l- 



ready close to the threshold value ([Acheson 1 119791 ). 



Holzwarth fc Schiissler ( 2003al lb) have studied tidal ef- 
fects on the stability and evolution of magnetic flux 
tubes in stellar interiors in the case of close binary 
systems. The maximum tidal perturbation occurs 
for flux tube s store d in the subphotospheric layers. 
Brandenbur j ( 20051 ) has recently discussed the pos- 
sibility that the solar dynamo has a source region 
located in the subphotospheric layers of the convec- 
tion z one, although other authors, e.g., iDikpati et al 



to operate there. Nevertheless, the possibility of a flux 
storage not far from the surface cannot be excluded in 
a magnetically active star. 

The main difficulty with this model is that tidal ef- 
fects are virtually negligible for CoRoT-4 and CoroT-6 
since the ratio of the semimajor axis of the orbit to 
the stellar radius a/R ^ 17—19 in those systems and 
the tidal deformation scales as {a/R)~^ . Moreover, two 
tidal bulges are expected, separated by 180° in longi- 
tude, but only one active longitude is associated with 
the planet in those stars. Therefore, it is unlikely that 
a tidally-based mechanism can account for the observa- 
tions. 

An a lterna tive mechanism has been suggested by 
Lanza I (I2OO8I ) in the appendix to his paper. It as- 
sumes that the stellar dynamo is somehow perturbed 
by the planet. Two main effects are responsible for 
the magnetic field generation and modulation in late- 
type stars: differential rotation and the a effect, the 
latter related to the mean helicity of the turbulent mo- 
tions. Differential rotation stresses poloidal fields di- 
rected along the meridians and produces a toroidal field 
component, i.e., parallel to the equator, that gives rise 
to the sunspots after emerging as discrete flux tubes. 
The a effect is responsible for the re-generation of 
the poloidal field from the toroidal field allowing the 
dynamo to operate. As discussed i n, e.g., Gh. 9 of 



Brandenburg fc Subramanian (|2005l) . also the mean 



current helicity associated with small-scale magnetic 
flelds plays a role in generating the a effect. The ex- 
pression for the a effect that takes into account the 
current helicity is: 



a oc {{v • V X t/) — (j • 6))t, 



(2) 



where v is the turbulent velocity field, h the turbulent 
magnetic field, j = V xb the turbulent current density, 
the angular brackets denote an ensemble average, and r 
is a characteristic correlation time between the turbu- 
lent fields. The mean current helicity is related to the 
evolution of the total helicity of the magnetic field of 
the star by the law of helicity conservation (s ee Ghs. 3 
and 9 of iBrandenburg fc Subramanian 1120051 ). Specifi- 
cally, the variation of the magnetic helicity of the star 
is given by: 



dt 



A-B dV 



-2t] j B dV, 
'v 



(3) 



(|2002l ). pointed out several difficulties for the dynamo 



where A is the vector potential of the magnetic field 
B = V X A, V the volume of the star, r] the magnetic 
diffusivity of the plasma assumed to be uniform, and 
J = V X B the current density. For an ideal plasma 
77 = and the integral in the l.h.s., i.e., the magnetic 
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Fig. 5 A sketch of the three main phases of the conjectured star-planet interaction in the outer stellar corona. In the 
first stage (left panel), the magnetic reconnection connects a tall coronal loop with the planetary magnetosphere in a region 
of weak field intensity {B ~ 0.01 G); in the second stage (middle panel), the interconnecting loop is stressed and bent by 
the orbital motion of the planet that is not synchronized with the stellar rotation; finally, in the third stage (right panel), 
the energy stored in the interconnecting loop is released when it reconnects with low magnetic loops close to the star in a 
region of relatively strong (B ~ 1 G) field. 



helicity, is conserved. In the presence of dissipation, 
the rate of change of the magnetic helicity depends on 
the difFusivity and the integral in the r.h.s., i.e., the 
current heUcity. If we assume that the timescale for the 
variation of the magnetic helicity is much longer than 
the timescale for the variations of the currents, i.e., the 
dynamo operates in an almost stationary regime, then: 

/ J-BdV^O. (4) 
Jv 

Considering the current density and the magnetic field 
as the sum of their average and fluctuating components, 
viz. J = {J) +j and B = 0) + 6, Eq. ^ and the 
Reynolds rules give: 

/ {J) ■ {B) dV ^ - f {j-b)dV. (5) 
Jv Jv 

Equation ([5]) implies that the loss of the mean cur- 
rent helicity occurring at the reconnection sites in the 
corona, that changes the integral in the l.h.s., produces 
a variation of the small-scale helicity on the r.h.s. This 
in turn perturbs the a effect. Specifically, the reconnec- 
tion events triggered by the planet in the corona (see 
Sect. 14. 2p may induce a non-axisymmetric component 
of the a effect modulated in phase with the passage of 
the planet across the large scale loops. This phase- 
dependent a effect can produce a non-axisymmetric 
magnetic field that emerges at the surface where the 



intensity of the unperturbed magnetic field is already 
close to the threshold for the onset of the buoyancy 
instability. 

In conclusion, the perturbation of the a effect at the 
longitude where magnetic reconnection is more effec- 
tive, owing to the interaction with the planet, leads to 
a modulation of the field emergence with the orbital pe- 
riod of the planet or its synodic period. The timescale 
for such a modulation depends on the intensity of the 
background field and its amplification timescale that 
are ruled by the unperturbed stellar dynamo. In the 
case of CoRoT-2, we conjecture that two interconnect- 
ing loops are responsible for the observed modulation. 
The one with the footpoints in the latitude range oc- 
culted by the planet is located in a zone of higher dy- 
namo efficiency leading to a typical timescale of 6 syn- 
odic periods for the magnetic flux modulation. On the 
other hand, the one located at higher latitudes, where 
the dynamo is less efficient, requires 10 synodic periods 
between successive emergence episodes because the dy- 
namo takes longer to bring the field close to the thresh- 
old. 



5 Discussion and conclusions 

We have reviewed the current evidence of star-planet 
magnetic interaction in stellar atmospheres focussing 
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on the photospheres. Unfortunately, the current re- 
suhs neither point definitely to a single type of phe- 
nomena nor can be explained by a single model, so fur- 
ther dedicated observations are needed to assess the 
reality and clarify the properties of the related phe- 
nomena. Nevertheless, one common aspect of the dif- 
ferent phenomenologies appears to be the presence of 
starspots whose formation is triggered by the passage 
of the planet over stellar active longitudes, and that are 
often localized close to the subplanetary longitude. 

Space-borne long-term photometry such as that 
made possible by MOST, CoRoT, Kepler (and later 
PLATO) provides unique data sets to investigate the 
presence of photospheric starspots related with a hot 
Jupiter, as we discussed in the cases of CoRoT-2, 
CoRoT-4, CoRoT-6, and Kepler targets. Ground-based 
searches should concentrate on the monitoring of chro- 
mospheric proxies, such as Ca II H & K line core emis- 
sion, to obtain a sufficient statistics on the properties 
of chromospheric hot spots. Those observations pro- 
vide basic constraints to the magnetic interaction in the 
outer atmospheric domains that may induce perturba- 
tions of the stellar dynamo leading to the formation of 
starspots associated with the planet (cf. Sect. 14. 3p . 

Theoretical models should clarify the mechanisms for 
the interaction and this requires a map of the coronal 
field structure. As in the case of the Sun, it is possible to 
extrapolate the field components measured at the pho- 
tosphere. Since the high-order multipoles of the photo- 
spheric field decay rapidly with the distance, only the 
low-order components, say dipole and quadrupole, are 
relevant at the distance of the planet. This means that 
the present photospheric field reconstructions based on 
spectropolarimetric techniques are adequate to extrap- 
olate the field at the distance of the planetary magneto - 



sphere (e.g., iMoutou et al. IboOTt [Pares et al. Iboiol) . 



The simultaneous observations of a chromospheric hot 
spot and the photospheric field to extrapolate the coro- 
nal configuration, will provide unique inforrn ation to 
understand the interaction (see Lanza 20091 for de- 
tails). 

The possibility that the planet somehow triggers the 
emergence of new magnetic flux from the interior, ap- 
pearing as starspots or faculae in the photosphere, can- 
not be excluded and seems to be supported by the ob- 
servations of the CoRoT and Kepler targets discussed 
above. If confirmed in other systems, this will open a 
new view on stellar dynamos, planetary magnetic fields 
and star-planet interaction making this area of research 
of the highest interest for the advancement of stellar 
and planetary astrophysics. 
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